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In this contribution two recent analyses for the extraction of the charm quark 
mass are discussed. Although they rely on completely different experimental 
and theoretical input the two methods provide the same final results for the 
charm quark mass and have an uncertainty of about 1%. 
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1. Introduction 

There has been an enormous progress in the determination of the quark 
masses in the recent years due to improved experimental results, many high- 
order calculations in perturbative QCD and precise lattice simulations.^ In 
this contribution we describe two recent analyses which lead to the most 
precise results for the MS charm quark mass. 

The first method^"^ is based on four-loop perturbative calculations for 
the moments of the vector correlator which are combined with moments ex- 
tracted from precise experimental input for the total hadronic cross section 
in electron positron collisions. 

Also the second method^ relies on four-loop calculations, however, for 
the pseudo-scalar rather than for the vector current correlator. It is com- 
bined with data obtained from simulations on the lattice with dynamical 
charm quarks. The latter are tuned such that the mass splitting between 
the T' and T and the meson masses m^, 2m^ — m^, m^^ and mx are cor- 
rectly reproduced. Thus the underlying experimental data are completely 
different from the first approach. 



September 11, 2008 8:1 



WSPC - Proceedings Trim Size: 9in x 6in caqcd 



2. R{s) and perturbative QCD 

The basic object for the first method is the total hadronic cross section in 
e+e~ annihilation. Normalized to the production cross section of a muon 
pair it defines the quantity 

, cr(e+e~ hadrons) 
R{s) = , (1) 

where apt = 47ra^/(3s). 

A compilation of the experimental data contributing to R{s) in the 
charm region can be found in Fig. 1. For our analysis it is of particular 
importance to have precise values for the electronic widths of the narrow 
resonances J/^f and which have been measured by various experiments.-'^ 
Furthermore, we rely on the excellent data provided by the BES collabo- 
ration^'^ in the region between 3.73 GeV (which is the onset of D meson 
production) and about 5 GeV which marks the end point of the strong 
variations of R{s). Above 5 GeV R{s) is basically flat and can be described 
very well within perturbative QCD taking into account charm quark effects. 
Thus in this region we use rhad,^ a f ortran program containing all state- 
of-the-art radiative corrections to R{s) since between 5 GeV and 7 GeV no 
reliable data is available. 

Since we are interested in the extraction of the charm quark mass we 
have to consider the part of R{s) which corresponds to the production 
of charm quarks, usually denoted by Rc{s). Rc{s) is used to compute the 
so-called experimental moments through 

^r^/^w. (2) 

It is clear that in order to perform the integration in Eq. (2) one has to 
subtract the contributions from the three light quarks. This has to be done 
in a careful manner which is described in detail in Ref.^ 
The theoretical counterpart to Eq. (2) is given by 

where the Cn are obtained from the Taylor coefficients of the photon po- 
larization function for small external momentum. 

Low moments are perturbative and have long been known through three- 
loop order ^^"^^ (see Ref.^^'^^ for moments up to n = 30). More recently also 
the four-loop contribution for n = 1^^'^'' and n = 2 could be evaluated^^ 
(see also Ref.^^). 
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Fig. 1. R{s) around the charm threshold region. The solid line corresponds to the 
theoretical prediction. The uncertainties, which are indicated by the dashed curves, are 
obtained from the variation of the input parameters and of /i. The inner and outer error 
bars give the statistical and systematical uncertainty, respectively. Next to the data from 
BES"''^ we also show the results form MD-1^ and CLEO.^" The narrow resonances are 
indicated by dashed lines. 



Ill the perturbative calculation we renormalize the charm quark mass 
in the MS scheme. This enables us to extract directly the corresponding 
short-distance quantity avoiding the detour to the pole mass and the cor- 
responding intrinsic uncertainty. 

The results obtained for the charm quark mass from equating the exper- 
imental and theoretical moments are collected in Tab. 1. In order to obtain 
these numbers we set the renormalization scale to = 3 GeV and extract 
as a consequence mc(3 GeV). The uncertainties are due to the experimental 
moments, 6as{Mz) — ±0.002, the variation of /i between 2 GeV and 4 GeV 
and the non-perturbative giuon condensate. 

Ill contrast to the corresponding table in Ref.'* we included in Tab. 1 
the new four- loop results from Ref.^^ for n = 2. This leads to a shift in the 
central value from 0.979 GeV to 0.976 GeV. Furthermore the uncertainty 
of 6 MeV which was due to the absence of the four-loop result is removed. 

The results in Tab. 1 show an impressive consistency when going from 
n = 1 to n = 4 although the relative weight form the various energy regions 
contributing to is completely different: whereas for n = 1 the region 

for y/s > 5 GeV amounts to about 50% of the resonance contribution it is 
less than 4% for n = 3. Also the decomposition of the uncertainty changes 
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Table 1. Results for mc(3 GeV) in GeV. The errors are from 
experiment, as, variation of fi and the gluon condensate. The 
error from the yet unknown four-loop term is kept separate. 



n 


■mc(3 GeV) 


exp 


as 




np 


total 




1 


0.986 


0.009 


0.009 


0.002 


0.001 


0.013 




2 


0.976 


0.006 


0.014 


0.005 


0.000 


0.016 




3 


0.982 


0.005 


0.014 


0.007 


0.002 


0.016 


0.010 


4 


1.012 


0.003 


0.008 


0.030 


0.007 


0.032 


0.016 



1.3 r 

1.25 - 




n 

Fig. 2. mc(3 GeV) for n = 1,2,3 and 4. For each value of n the results from left to 
right correspond the inclusion of the one-, two-, three- and four-loop terms in the theory 
moments. 

substantially as can be seen in Tab. 1. Whereas for n = 1 the contribution 
from the fi variation is negligible it exceeds the experimental uncertainty 
for n — 3. 

In Fig. 2 we show for the first four moments the result for rric (3 GeV) as 
a function of the loop order used for A4l^. One observes a nice convergence 
for each n. Furthermore, the consistency among the three- and in particular 
the four-loop results is clearly visible from this plot. 

As final result of the analysis described in this Section we quote the 
value given in Ref.* which reads 

mc(3 GeV) = 0.986(13) GeV . (4) 
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Table 2. Results for mc(3 GeV) in GeV. Both the to- 
tal uncertainties are shown and the splitting into contri- 
butions from the lattice simulation, as, missing higher or- 
der corrections and the non-perturbative gluon condensate. 



n 


mc(3 GcV) 


lattice 


as 


h.o. 


np 


total 


2 


0.986 


0.008 


0.003 


0.004 


0.003 


0.010 


3 


0.986 


0.009 


0.004 


0.003 


0.000 


0.011 



3. Lattice gauge theory and perturbative QCD 

In the recent years there has been a tremendous progress in developing 
precise QCD simulations on the lattice. In particular, it has been possible 
to simulate relativistic charm quarks using the so-called Highly Improved 
Staggered Quark (HISQ) discretization of the quark action. ^"'^^ In Ref.^ 
this has been used to evaluate moments of the pseudo-scalar correlator 
with an uncertainty below 1%. The moments from the lattice calculation 
arc equated with the ones computed within perturbative QCD. In Rcf.-'^^ 
the second non-trivial moment could be evaluated with the help the axial 
Ward identity from the first moment of the longitudinal part of the axial- 
vector current. Very recently this trick could be extended in order to arrive 
at the third moment for the pseudo-scalar current. 

Tab. 2 summarizes the results obtained for mc(3 GeV) (for n = 2 and 3)'^ 
together with the corresponding uncertainties from the lattice, ag , missing 
higher order perturbative corrections and the gluon condensate.*' 

Like in the previous section we find also here an excellent agreement in 
the central values which leads us to the final result 

mc(3 GeV) = 0.986(10) GeV. (5) 

Let us mention that the dimonsionless first moment can be used to 
extract a value for the strong coupling. We can furthermore consider ratios 
of moments in order to get rid of the overall dependence on nic and again 
extract a^. In Ref."'' this has been done for the ratio of the second to the 
third moment which is known to four-loop order within perturbative QCD. 
The two determinations lead to 

(3 GeV) =0.251(6), (6) 



''Note that for n = 1 no charm quark mass can be determined since, in contrast to the 
vector correlator, the corresponding moment is dimensionless. 

''For the presentation in this Section the notation of Ref.® for the numeration of the 
moments has been translated to the one of Ref.* 
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which corresponds to'^ 



a' 



;(5)(M2) = 0.1174(12). 



(7) 



This value agree well with the particle data group result and other recent 
determinations (see, e.g., Refs.^**'^^). 

4. Summary 

In this contribution we have presented the two to date most precise de- 
terminations of the charm quark mass. Let us stress once again that, al- 
though in both cases moments of current correlators are considered, the two 
methods rely on completely different experimental input and on different 
theory calculations. Whereas in one case perturbative QCD is compared 
with experimental data for i?(s), in the second case high precision lattice 
simulations with dynamical charm quarks are crucial ingredients. It is quite 
impressive that the final results as given in Eqs. (4) and (5) coincide both 
in the central value and the uncertainty. 

In Fig. 3 we compare the results of Section 2 and Section 3 with various 
other recent determinations. One observes a good agreement, however, our 
results are by far the most precise ones, as can be secin by the grey band. 

Up to this point we have presented results for the MS charm quark 
mass evaluated at the scale = 3 GeV. In general, the comparison of 
results from various analyses are performed for the scale- invariant mass, 
Tnc{mc) (sec, e.g., Rcf.-'^). Note, however, that the scale /i = rric is quite 
low and the numerical value of is relatively big. Thus, it would be more 
appropriate to perform the comparison at a higher scale like fj, = 3 GeV. 
Let us nevertheless present the scale-invariant charm quark mass. From 
0.986(10) GeV one obtains 



The method described in Section 2 can also be used to extract the 



After including the new four- loop results from Ref.^^ the result of Eq. (9) 
becomes 



mc(mc) = 1.268(9) GcV . 



(8) 



bottom quark mass. The analysis of Rcf.'* leads to 
nibimb) = 4.164(25) GcV. 



(9) 



mb{mb) = 4.162(19) GeV 
which has a significantly reduced uncertainty. 



(10) 



"^The calculation of the running and decoupling is easily done with the help of RunDec. 
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Fig. 3. Comparison of recent determinations of mc(3 GeV). 
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